Quasiparticle Localization in Disordered d-Wave Superconductors 
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An extensive numerical study is reported on disorder effect in two-dimensional d-wave supercon- 
ductors with random impurities in the unitary limit. It is found that a sharp resonant peak shows 
up in the density of states at zero energy and correspondingly the finite-size spin conductance is 
strongly enhanced which results in a non-universal feature in one-parameter scaling. However, all 
quasiparticle states remain localized, indicating that the resonant density peak alone is not sufficient 
to induce delocalization. In the weak disorder limit, the localization length is so long that the spin 
conductance at small sample size is close to the universal value predicted by Lee (Phys. Rev. Lett. 
71, 1887 (1993)). 
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Since the discovery of the ci-wave pairing symmetry 
in high-T c cuprates, there has been increased interest in 
low-energy quasiparticle properties in unconventional su- 
perconductors. In such a d-wave superconductor, quasi- 
particles are gapless along the four nodal directions on 
the essentially cylindrical Fermi surface, in contrast to 
the conventional s-wave superconductors, where quasi- 
particles are gapped everywhere. The issues of how the 
disorder affects the low energy quasiparticle excitations 
and whether these quasiparticles are localized remain 
unresolved. Some perturbative self-consistent T-matrix 
(SCTM) calculations @-|] and a non-perturbative one M 
predicted a nonzero constant density of states (DOS) in 
low energy region in the presence of weak disorder. How- 
ever, most of non-perturbative calculations Jl(|[llJ and 
one numerical study 1 12 showed that the DOS at zero en- 
ergy vanishes. With this constant or even vanishing DOS, 
some groups [||,[l3| suggested that all quasiparticle states 
are localized. On the other hand, it has been shown |l4| ] 
that a single unitary scattering impurity produces a zero- 
energy quasiparticle resonant state while the long-range 
overlap between these impurity states |l5| may lead to 
extended quasiparticle band near zero energy. More re- 
cently, a singularity in the DOS at zero energy was ob- 
tained by non-perturbative T-matrix method fl6| for the 
random distributed unitary impurities. It is noteworthy 
that in one dimension, there is a direct relation between 
the localization length and the DOS p7| ; thus a singu- 
larity in zero-energy DOS signals the delocalization in 
the system. However, in two dimensions, this theorem 
does not hold generally |i"rj| ]. For example, in two dimen- 
sional integer quantum Hall systems, it has been shown 
that the delocalization property at the quantum critical 
point is not changed by the changing of the DOS due to 
strong electron-electron interaction |fL8| . Therefore, the 
localization of quasiparticles in a ci-wave superconductor 
in the presence of nonmagnetic unitary impurities is still 
an open question. 

In this Letter, we numerically examine the disorder 
effect in d-wave superconductors with nonmagnetic im- 



purities in the unitary limit. The quasiparticle DOS is 
calculated by exact diagonalization and the spin conduc- 
tance is computed by the transfer matrix method. It 
is found that, depending on the particle- hole symme- 
try of the Hamiltonian, a sharp DOS peak can occur 
at zero energy and correspondingly the spin conductance 
is strongly enhanced at finite sample size. However, us- 
ing one parameter scaling analysis we show that all the 
quasiparticle states are always localized regardless of the 
existence of the zero-energy peak in the DOS. In weak 
disorder limit, the localization length is so long that the 
spin conductance at small sample size remains close to 
the universal value 2£,o/a (£o is the coherence length of 
the superconductor and a is the lattice constant) in agree- 
ment with the theoretical prediction by Lee [Q. A non- 
universal feature in one-parameter scaling of conductance 
related to the resonant DOS peak is also discussed. 

We begin with a lattice Hamiltonian for the d-wave 
superconductor |l9| 

where (ij) refers to two nearest neighboring sites with 
the hopping integral taken as the unit, fi is the chemical 
potential and Ui is the impurity potential. We mainly 
consider the unitary limit where Ui takes a nonzero 
value Uq only at a fraction n.; of the sites which are 
randomly distributed in space. The d-wave symmetry 
is imposed by choosing order parameters: A^i±$ — 
—Ai^±y — Ad, which yields the excitation spectrum 

-2(cosfc x + cosky) — n and 
cosfcy). Therefore, gapless quasipar- 
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ticle states exist along the direction \k x \ — \k y \ in the mo- 
mentum space. Unless otherwise stated, fi = 0, A^ = 0.1, 
and Uo = 100 are taken throughout the work. In the 
presence of a single impurity in the unitary scattering 
limit, earlier study has shown ]20[ ] that the order param- 
eter is strongly suppressed near the impurity site on a 
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scale of a few lattice constant. To take into account this 
effect, the order parameters on bonds connecting with 
the strong impurity sites are taken as zero. By exactly 
diagonalizing the Hamiltonian (1), one can calculate the 
quasiparticle DOS, which is defined as ]l^]2~lj ] 

P{E)^^- L Y J 5{E-E n ), (2) 

n 

where Nl — L x L with L the linear dimension of the sys- 
tem in units of lattice constant (a = 1). In Fig. |l](a), the 
DOS is plotted as a function of quasiparticle energy E 
with impurity density rii = 0.1 (/i = and = 0.1) at 
Nl = 90 x 90. In the calculation the periodic boundary 
condition is used and it has been checked that the results 
do not depend on the boundary condition for large L con- 
sidered here. As shown in Fig. 0(a), we find that a sharp 
zero-energy peak shows up in the DOS, which can be fit- 
ted by the analytical form c ni/2\E\(\n 2 \E/E g \ + 7r 2 /4) 
from Ref. fl6|| with c = 0.66 (the superconducting gap 
Eg ~ 4 x A d ) as shown in the inset of Fig. |](a). The 
fitting breaks down at an energy scale close to 1 /Nl [§2| . 
The strength of the zero-energy peak is reduced when n, 
is changed to 0.04 as shown in E < part of Fig. [|(b). 
The overall shape of the peak is sample size independent 
(from Nl = 25 x 25 to 120 x 120 as well as a strip system 
30 x 300) with the DOS value at the peak position in- 
creasing with Nl very slowly. Note that in the presence 
of disorder, the order parameter Ay is in principle sub- 
ject to the self-consistency condition: Ay = —gij{cj±Cii) 
where <jy = go is the attracting interaction for d-wave 
pairing. We have also calculated the DOS by diagonal- 
izing the Hamiltonian Eq. (jlj) self-consistently for each 
disorder configuration with Nl = 25 x 25. The obtained 
DOS is also shown in Fig. 1(b) (E > part), which is av- 
eraged over 50 impurity configurations and 8x8 wavevec- 
tors in the supercell Brillouin zone. As can be seen, the 
DOS value at E g is reduced due to the suppression of Ay 
around each impurity site. However, all other features re- 
main essentially unchanged (compared to the E < part 
of Fig. 1(b)). 

The presence of zero-energy peak in the DOS crucially 
depends on the symmetry of the Hamiltonian. Since 
a repulsive or attractive impurity center with infinite 
strength (i.e., unitary limit) is equivalent to the exclusion 
of a lattice site, both the local and global particle-hole 
symmetry remain if the chemical potential a = 0, which 
then produces a resonant peak at E = ]lq| . This fea- 
ture was not exhibited in earlier works because of cither 
the breaking of local particle-hole symmetry by the soft 
impurity scattering |l^,|ll|,|l3| or the breaking of the band 
particle- hole symmetry by considering /i ^ fOg] , which 
indicates the importance of the realization of disorder 
model. In Fig. |l|(c), the DOS is presented for random 
on-site disorders with Ui uniformly distributed between 
[-1,1] with a homogeneous order parameter at all sites. 



Due to the absence of the local particle-hole symmetry 
in this case, p{E) has a finite value at low energy region 
down to a mesoscopic scale E ~ 1/Nl, below which it 
shows a zero-energy dip, in agreement with those of the 
non-perturbative calculations |T(],[ll],[l3|] . We have also 
relaxed the unitary scattering limit by taking Uq = 10 
(comparable to the band width) or (and) broken the band 
particle- hole symmetry by taking ji = — 1, and found that 
the DOS at E — is always strongly suppressed, which 
is similar to the results obtained in Ref. p2| . There is a 
smooth crossover from the zero-energy peak to dip in the 
DOS with the varying of model parameters as long as the 
Hamiltonian is driven away from the perfect particle-hole 
symmetry. Therefore, our numerical result indicates that 
different realizations of disorder give rise to different pro- 
files of the DOS as the energy approaches to the Fermi 
level. 

In the absence of the zero-energy resonant peak in DOS 
due to the breaking of band particle-hole symmetry by 
fi ^ 0, Franz et al. |23j | have studied the similar prob- 
lem by examining the sensitivity of the wave function 
to the boundary conditions and by analyzing the finite- 
size dependence of inverse participation ratios, and pre- 
sented a strong evidence for the localization of low en- 
ergy quasiparticles. Here we are concerned with the ques- 
tion of quasiparticle localization or derealization in the 
disordered d-wave superconductor when the zero-energy 
peak appears in the DOS. We employ the transfer-matrix 
method to calculate the finite-size localization length and 
the longitudinal conductance. We consider the quasi-one- 
dimensional strip sample with the length L > 10 5 and 
width M. The quasiparticle wavefunction amplitudes in 
the ix-th and (ix + l)-th slices satisfy the following equa- 
tion: 




where (j>i X is a 2Af-component vector of the Bogoliubov 
amplitudes for quasiparticle states, and Ti X is a 4M x AM 
transfer matrix. The transfer matrix through the whole 
system, Pl = Y]_ ix= \ T ix , has a set of 2M pairs of Lya- 
punov exponents, which determine the inverse of length 
Xi (i = 1, 2, . . . , 2M). The usual orthonormalization pro- 
cedure is taken p4j in our calculation. Correspondingly, 
the longitudinal conductance g s extrapolated for square 
sample with width M is given by p5[ : 

2M 

g s {M,m) = ^cosh _2 Aj , (4) 
j=i 

where Aj = Xj/M. Note that g s corresponds to the spin 
conductance as the spin carried by quasiparticle is con- 
served As shown in Fig. ||(a), g s as a function of M 
monotonically decreases with the increase of M at E = 
for each selected impurity density from m = 0.01 to 0.16, 
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consistent with localization in the large M limit. All the 
data between M = 32 and M = 120 at different m can 
be collapsed onto a single curve: 



g s (M,m) = f 



M 



(5) 



as shown in Fig. g(b) in accordance with the one- 
parameter scaling law. Here is the thermodynamic 
localization length which only depends on ni as shown 
in the inset of Fig. ||(b), and it remains finite for all the 
disorder density rij, suggesting that all the states are lo- 
calized even in the unitary limit with the presence of the 
zero energy resonant peak. In addition, we display in 
Fig. ^| the conductance g s as a function of quasiparticle 
energy E at sample width M — 48 and 96 with m =0.1. 
It is found that g s is strongly enhanced as E — > in 
the region E ~ 0.01 corresponding to the width of the 
DOS peak while away from this region g s generally in- 
creases with the increase of E. For quasiparticle at low 
energy with E < 0.1, it has been found that g s always 
decreases with the increase of M and all the states are 
localized in large M limit. However, the scaling curve 
g s (M,m) = f(^p-) found at E away from the DOS 
peak (E > 0.01) is different from that for E = at the 
peak of the DOS, indicating the breaking-down of the 
universal one-parameter scaling law due to the presence 
of the resonant peak in the DOS. In addition, for all val- 
ues of Ad, <7s decreases with M, implying localization in 
all the parameter region. At fixed M, g s always decreases 
monotonically with the increasing Ad, in agreement with 
the general argument that quasiparticle states in a super- 
conducting phase are always more localized [Q than the 
corresponding normal state. Given the localization of the 
quasiparticle states, the spin conductance at the longest 
length scales must vanish. However, in the weak disorder 
(the impurity density rn < 4%), i.e., the Born limit, we 
found that the localization length is so long that at small 
sample size M ~ 32, the spin conductance g° follows the 
universal form Q g° = 2£o/a as the coherence length 
£o is changed from 1.6 to 6.4 by changing Ad between 
0.05 and 0.20. For the strong disorder, the conductance 
is generally smaller than the Born limit value due to the 
onset of localization effect. 

In conclusion, we have studied the quasiparticle states 
in 2D d-wave superconductors with randomly distributed 
strong impurities in the unitary scattering limit. As the 
particle-hole symmetry holds, a very sharp DOS peak is 
obtained at zero energy. Such a DOS peak enhances the 
finite-size conductance. However, using one-parameter 
scaling analysis we have shown that all the quasiparti- 
cle states at low energy are still localized and the local- 
ization effect is generally enhanced with the increase of 
superconducting order parameter. 

Note added: After the submission of this paper, we 
received a preprint from Atkinson et al. [p6[ where similar 



results for the DOS in the unitary limit of the symmetric 
band were obtained. 
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FIG. 1. (a) Density of states p(E) in a d-wave superconduc- 
tor (sample size 90 x 90) with a fraction m = 0.1 of the ran- 
domly distributed impurities in the unitary scattering limit 
Uo = 100. Inset of (a): The DOS peak is fitted using the 
analytical form con;/2|£|(ln 2 \E/A g \ + tt 2 /4) (E g = 4 x A d ) 
from Ref. (l(J with c = 0.66. (b) Comparison of p(E) (E < 
part) at Nl = 90 x 90 for rii = 0.04 with the self-consistent 
result (E > part) at Nl — 25 x 25 (averaged over 8x8 
supercells). (c) p(E) at Nl = 90 x 90 with on-site random 
disorders distributed between [— 1, 1]. 

FIG. 2. (a) Spin conductance g s as a function of the strip 
width M for different impurity density rii at E — 0. (b) 
Double logarithmic plot of g s as a scaling function of M/£ at 
E = for all the data with 0.01 < m < 0.275. Inset of (b): 
The scaling parameter £ as a function of rii at E — 0. The 
other parameters are p — 0, Ad = 0.1, and Uo ~ 100. 

FIG. 3. g s as a function of quasiparticle energy E at 
M — 48 and 96 for rii = 0.1. All the other parameters are the 
same as in Fig. 2. 
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